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■ ABSTRACT 

(N 

Context. 

Aims. We present a new method that uses luminosity or stellar mass functions combined with clustering measurements to select 
samples of galaxies at different redshifts likely to follow a progenitor-to-descendant relationship. As the method uses clustering 
information, we refer to galaxy samples selected this way as clustering-selected samples. We apply this method to infer the number 
of mergers during the evolution of MUSYC early-type galaxies (ETGs) from z ^ 1 to the present-day. 
(~| ■ Methods. The method consists in using clustering information to infer the typical dark-matter halo mass of the hosts of the selected 

progenitor galaxies. Using ACDM predictions, it is then possible to follow these haloes to a later time where the sample of descendants 
I I will be that with the clustering of these descendant haloes. 

O ■ Results. This technique shows that ETGs at a given redshift evolve into brighter galaxies at lower redshifts (considering rest-frame, 

passively evolved optical luminosities). This indicates that the stellar mass of these galaxies increases with time and that, in principle, 
C/3 ' a stellar mass selection at different redshifts does not provide samples of galaxies in a progenitor-descendant relationship. 

I Conclusions. The comparison between high redshift ETGs and their likely descendants at z = points to a higher number density for 

the progenitors by a factor 5.5 ± 4.0, implying the need for mergers to decrease their number density by today. Because the luminosity 
densities of progenitors and descendants are consistent, our results show no need for significant star-formation in ETGs since z = 1, 
^ which indicates that the needed mergers are dry, i.e. gas free. 

' Key words, galaxies: formation and evolution - statistics - cosmology: observations 

\o 

1. Introduction gies, such as the A-Cold Dark Matter (ACDM) model, have 

. been compared to observational data from different redshifts 

O The study of the population of early type galaxies (ETGs) at dif- wavelength ranges, via numerous statistical probes; in some 

^ ferent redshifts has been used extensively to test our knowledge ^^^^^ ^j^^ agreement appeal's to be excellent and in other cases 

T-H , of the galaxy formation process, in particular, the assembly of jj^^^^ ^^^^ important discrepancies. Some earUer failures 

^ . the stellar content of massive galaxies. One of the main advan- ^j^^ ^^^^^j^ j^^^^ ^^^^ consequences of adopting the 

■ ' tages in studying ETGs resides in their simple evolution, which ^ cosmological model (see for instance Heyl et al., 1995), 

allows to use passive evolution to infer their likely z = lu- ^jj^^^.^ jj^^ j^^,^ important astrophysical processes in 

H ■ minosities taking into account the aging of their stellar popula- ^j^^ formation recipe, such as for instance, galaxy merg- 

y . tions; this allows to use passively evolved luminosities (PEL) as ^ ^^^^ ^994) and feedback from AGN which 

good proxies for their stellar mass. Analyses of the evolution of ^j^^^ j^^j^^^^^ ^^1^^^ discrepancies in the number density of 

the stellar mass and PEL functions have found that high stellar ^^.^ j^^^^^ ^ B^^^^ 2006). More recent models 

mass (M, > lO^h-'Mo), passive galaxies do not show evolu- ^j^^^ reasonable agreement with galaxies in the local Universe 

tion m their comoving space density since z ~ 1 (Cimatti et (Kauffmann et al., 1999; Cole et al., 2000; Bower et al., 2006; 

al., 2002, 2004; McCarthy et al., 2004; Glazebrook et al. 2004; ^ucia et al., 2006; Lagos, Cora & Padilla, 2008; Khochfai- 

Daddi et al., 2005; Saracco et al., 2005; Bundy et al., 2006; ^ 0^^^.;,^^^^ 2008; Lagos, Padilla & Cora, 2009), but still strug- 
Perez-Gonzalez et al., 2008; Marchesini et al., 2009; Ferreras 



X 



gle to agree with a number of recent statistical measurements of 
et al., 2009). This result has been interpreted as evidence that , ^ 1 ^o 3 galaxies, such as the stellar mass function (Perez- 
the stellar content of such galaxies is already m place at high G^n^alez et al., 2008; Marchesini et al., 2009; but see Benson 
redshifts ruhng out the involvement of mergers (even dry i.e., ^ Cevereux, 2010), the galaxy luminosity function (Marchesini 
gas-free) since z ~ 1 . ^^^^ Dokkum, 2007), and SCUBA (Holland et al., 1999) num- 

The details of the evolution of ETGs are still in debate 



Galaxy formation models embedded in hierarchical cosmolo- 



ber counts of z ^ 3 galaxies detected in sub-millimetre bands, 
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so far only succesfully reproduced by the Baugh et al. (2005) 
model. A slightly better agreement is found for z ~ 1 galaxies as 
in the comparison between model and observed Luminous Red 
Galaxies by Almeida et al. (2008). 

Hopkins et al. (2008) presented an analysis of the evolution 
of the density of ETGs from a theoretical point of view. They 
show that part of the apparent mismatch between models and 
observations could originate from an inadequate mechanism to 
quench the star-formation in galaxies which have experienced 
major mergers. Due to this neither the Croton et al. (2006) nor 
the Bower et al. (2007) star-formation quenching approaches, 
via a dark-matter host-halo mass threshold or disc instabilities, 
respectively, are able to reproduce the z ~ 1-3 measurements 
of the number density of bright galaxies. Hopkins et al. suggest 
an intermediate solution may fit observational data. 

It could also be possible that the disagreement between mod- 
els and observations is produced by systematic biases from large 
redshift measurement errors or selection effects afl'ecting the ob- 
servational statistics at high redshifts, or simply due to the difl^er- 
ent and unknown nature of high redshift galaxies. Some results 
indicate that the effect from large redshift errors does not seem to 
be a major factor. For instance the analysis of high redshift clus- 
ters in the infrared by Mancone et al. (2010) shows an evolution 
of 1 .5 magnitudes between redshifts of 1 and 0, which is con- 
sistent with passive evolution, even when clusters allow a better 
handle on redshift estimate errors. This conclusion was also ob- 
tained by Muzzin et al. (2008) for the GCLASS cluster survey at 
z=l. On the other hand, Marchesini et al. (2010) showed that al- 
lowing for the existence of a previously unrecognized population 
of massive, old, and very dusty galaxies at z ~ 2.6 would result 
in good agreement between the observed abundance of massive 
galaxies at z = 3.5 and that predicted by semi-analytic models. 

One common method to study the evolution of ETGs uses 
PEL function (PELF) or stellar mass function measurements of 
ETGs to follow the evolution of the number density of PEL 
or stellar mass-selected galaxies. This method was applied to 
catalogues with either photometric redshifts (photo-zs) such as 
COMBO- 17 (Bell et al, 2004, B04) and the Subaru/XMM- 
Newton Deep Survey (SXDS, Yamada et al., 2005), or with spec- 
troscopic redshifts such as DEEP2 (Faber et al., 2007). These 
LFs have been analysed by Cimatti, Daddi & Renzini (2006, 
CDR) who find that the ETG population is consistent with no 
evolution in their number density since z ~ 0.8 (see also Banerji 
et al., 2010). We will apply this method to the PELF measure- 
ments obtained from the Multi- wavelength Survey by Yale-Chile 
(MUSYC, Gawiser et al., 2006) by Chrisdein et al. (2009, C09). 
Our approach takes into account the possible systematic effect 
imprinted on the luminosity function measurements by errors in 
the photo-zs. This new measurement can help to study the pos- 
sible effects from cosmic variance on this particular statistics, 
essential to compare with results from galaxy formation models. 

The number density evolution of samples of ETGs of equal 
stellar mass cannot, in principle, be used to infer a sequence 
of mergers - or the lack thereof - since they do not necessar- 
ily constitute causally connected populations. This is simply due 
to the possibility that as time passes, the stellar mass of a galaxy, 
even if it is an early-type, can change. For instance, Robaina 
et al. (2010) use stellar mass-selected samples as the same (sta- 
tistically) evolving population, and by measuring the number of 
close pairs from correlation function estimates, they infer a small 
number of 0.7 mergers per galaxy since z = 1.2 (see also Patton 
et al., 2000, Le Fevre et al., 2000, Lin et al., 2004, Kartaltepe 
et al., 2007); these estimates could be affected by biases due to 



the possible lack of a descendant relationship between samples 
at different redshifts. 

The main focus of this paper is to devise and apply a new 
method of selection of galaxy samples so that these follow a 
progenitor-descendant relationship, providing a new, detailed 
test for galaxy formation models. In this first application of the 
method we combine the LF measurements with recent results 
on the clustering of ETGs from the MUSYC survey (Padilla et 
al., 2010, PIO), which indicate that the descendants of ETGs are 
not simply ETGs of similar stellar masses (stellar mass-selected 
samples) or evolved luminosities at lower redshifts. 

Our approach will broaden the range of phenomena that can 
be constrained by combining LF and clustering measurements 
since, for instance, Lee et al. (2009) already used this com- 
bination to estimate the star-formation duty cycle and how its 
efficiency depends on halo mass for star-forming galaxies at 
z - 4 to 6. Using the halo model (see for example Jing, Mo & 
Borner, 1998, Peacock & Smith, 2000, Cooray & Sheth, 2002), 
and also combining LF and clustering measurements of red- 
sequence galaxies in the Bootes field. Brown et al. (2008, see 
also Zheng, Coil & Zehavi, 2005, White et al., 2007, Wake et 
al., 2008, Matsuoka et al., 2010) were able to follow the growth 
of the stellar mass of central and satellite galaxies between red- 
shifts z = 1 and z = for dark-matter haloes of a given mass. In 
particular, our work will be complementary to this effort, since 
we will be able to follow direct descendants while not restricting 
our analysis to stellar mass-selected samples. 

In the following section we describe the LF measurements 
we use in our analysis; in Section 3 we follow the evolution 
of the number density of MUSYC ETG galaxies using PEL- 
selected samples. In Section 4 we present the method of cluster- 
ing selection that allows the construction of samples that follow 
a progenitor-to-descendant relationship, and use them to study 
the role of mergers in their evolution. We conclude in Section 5. 

2. ETG passively evolved luminosity functions 

We will analyse the evolution of ETGs using their measured LF 
by C09 from MUSYC. This survey comprises over 1.2 square 
degrees of sky imaged to 5cr AB depths of U,B,V,R = 26, 
/ = 25, z = 24 and J, K(AB) = 22.5, with extensive follow- 
up spectroscopy. The source detection is done on the combined 
BVR image down to a magnitude of 27 (see Gawiser et al., 2006, 
for further details). 

C09 proposed and applied a new technique, the photomet- 
ric maximum hkehhood (PML) method, to a subset of MUSYC 
comprising the Extended Chandra Deep Field South (ECDF-S), 
covering approximately 0.25 sq. degrees on the sky. The PML al- 
gorithm was used to measure the underlying luminosity function 
of galaxy populations characterised by different spectral types. 
The latter are parameterised with a set of SED templates from 
Coleman, Wu & Weedman (1980, CWW) or fixed superposi- 
tions of two CWW templates, extended into the UV regime us- 
ing Bruzual & Chariot (1993) models. For the present analysis, 
we will only use the LF corresponding to the two earliest type 
templates in this set, which correspond to an elliptical galaxy 
and to a E-(-20%Sbc mix. PIO demonstrate that the sample of 
ETGs selected this way is comparable to a selection of the red- 
sequence at each redshift (as adopted in e.g. B04, CDR, Brown 
et al., 2008). They also indicate that these ETGs are compatible 
with a sample resulting from &K < 22.5 selection. 

The main advantage of the PML method is that it has been 
designed to deal with photometric fluxes alone; it does not re- 
quire estimates of distances or intrinsic galaxy luminosities, and 
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Fig. 1. Luminosity functions of ETGs at different redshifts passively evolved to z = (with redshift increasing from the left to the 
right panels as indicated in the key). The top and bottom rows show estimates in the r- and B-band, respectively, and for comparison 
the black solid triangles show the z - 0.165 SDSS ETG LF from Benson et al. (2007) in all the top panels, and from B04 in all 
the bottom panels. Red circles show the results from COMBO-17 (B04) and blue open triangles those from DEEP2 (Faber et al., 
2007). The thick light-blue solid and dotted lines show results from MUSYC, obtained via the PML and a photo-z based maximum 
likelihood methods by C09, respectively. The amount of evolution applied to the B and r-band luminosities is indicated in each 
panel. The shaded areas delimit the bright and faint ETG populations (dark and light grey, respectively). The vertical dashed lines 
(only visible in the subpanels corresponding to z > 0.74) mark the MUSYC completeness limit. 



it does not need to assume a distribution function for redshift 
errors of any form. This is achieved by comparing a trial LF 
to the observed galaxy sample in a parameter space consist- 
ing of observed fluxes, whose error function is well understood. 
Therefore, the PML method explicitly avoids a second convolu- 
tion with the photo-z errors, being in principle free of this possi- 
ble source of systematic effects. 

For comparison purposes, C09 applied a maximum likeli- 
hood method to provide estimates of photo-zs with competi- 
tive accuracies with respect to other photo-z measurement al- 
gorithms such as BPZ (Benitez, 2000). In the case of MUSYC, 
the resulting redshift uncertainty is Az ^ 0.065 x (1 + z) (nor- 
maUsed median absolute deviation, see Hoaglin, Mostelen & 
Tukey 1983). C09 also used these estimates of photo-zs to 
measure the LF of ETGs using a classic maximum likelihood 
method. We will compare the results on ETG evolution using 
both C09 LF estimates. 

We use the rest-frame r- and B-band LF estimates from 
MUSYC ETGs to allow a direct comparison to LF measure- 
ments in both bands. All the ETG LFs are passively evolved 



down to z = by applying an empirical passive evolution recipe 
whereby the evolved luminosity can be found via 

Mb(z = 0) = Mb(z)-h 1.1 5z, (1) 

following results from van Dokkum & Stanford (2003), Treu et 
al. (2005), and di Serego AHghieri et al. (2005). PIO showed that 
this passive evolution recipe is well followed by a [Fe/H] - 0.3 
single stellar population (SSP) which at z = 1 is 3.5Gyr old, and 
used this SSP to work out the equivalent recipe for the r-band, 
which is fit by 

M,(z = 0) = Mriz) + 0.98z. (2) 

This SSP is characterised by a colour which varies slowly with 
redshift, 

B-r 1.32-0.17XZ. (3) 

We apply these evolution corrections to B- and r-band LF 
measurements from relatively high redshift samples and com- 
pare them to results at z - 0.165 from the Sloan Digital Sky 
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Survey (SDSS; York et al., 2000) for ETGs by B04 in the B- 
band (selected using the red-sequence), and in the r-band by 
Benson et al. (2007, Be07), who selected ETGs as sources with 
a dominant bulge component (an alternative ETG LF measure- 
ment from SDSS is provided by Ball et al., 2006). We find that 
the two estimates are perfectly consistent with one another when 
applying Eq.[3] When correcting the luminosities by passive evo- 
lution and using this quantity to select samples of ETGs (with 
little or no star formation), we effectively mimic a stellar mass 
selection of galaxies (CDR). 

The comparison between the ETG LFs at different redshifts 
is shown in Figure [J where the Be07 and B04 SDSS ETG LFs 
are plotted in all the top and bottom panels, respectively. The 
MUSYC PML and classic maximum likelihood measurements 
of the ETG LFs for different redshift ranges (indicated in the 
panels) are shown for the r- and B-band (top and bottom pan- 
els, respectively). The bottom panels also show the LF estimates 
from COMBO-17 and DEEP2 for the same redshifts. As can be 
seen, the bright ETG population from COMBO 17 matches the 
ETG SDSS result for all the redshift slices shown in the fig- 
ure. Only the high redshift DEEP2 data shows a lower comov- 
ing number density than the SDSS LF at this luminosity range. 
Therefore, from these results one would conclude that there is no 
evolution in the number density of evolved ETGs from z ~ 0.85 
onwards. Q The MUSYC results (light-blue lines) for the bright 
population, as obtained from the PML algorithm or even from 
the classic maximum likelihood method, on the other hand, show 
indications of some evolution in the number density of bright 
ETGs with redshift in both bands, even at z < 0.45. This is 
more clearly seen in the r-band, than in the B-band when com- 
paring to the z = SDSS LF since, in the latter case, the SDSS 
ETG LF characteristic luminosity appears to be slightly fainter 
than for the low redshift MUSYC ETG, even when the B- and r- 
band SDSS LFs are entirely consistent when taking into account 
Eq.[3] This could indicate that the MUSYC ETG samples at low 
redshift are slightly redder than the SDSS ETGs. We will leave 
the discussion of the faint population to the following section. 

3. Evolved ETG luminosity selection of samples in 
MUSYC 

We study the number density evolution of PEL-selected samples 
from MUSYC, and compare it to previous results from DEEP2, 
COMBO-17 and SXDS. 

This approach has been applied in several ocasions in the lit- 
erature, and we apply it in this section to MUSYC ETG galaxies. 
This procedure can only be used to follow the evolution of the 
number density of ETGs of a particular PEL (or stellar mass). 
Due to their simplicity, selections of this type have been assumed 
to provide samples of the same evolving population at differ- 
ent redshifts (e.g. CDR, Robaina et al., 2010, Marchesini et al., 
2009). However, they cannot be used to study the frequency of 
mergers in principle, since the stellar mass or evolved luminosity 
of an early-type galaxy can evolve with time, i.e., there is no pre- 
cise information on whether the two populations are connected 
in a progenitor-descendant relationship (PIO). 

Notice, however, that this approach still allows a comparison 
with models, and in particular it has been used in galaxy forma- 



Notice that this does not prove there are no mergers since it would 
still be possible that galaxies of a given luminosity merge to form 
brighter galaxies while maintaining a close to constant luminosity func- 
tion (see Section 2) with galaxies in this bin leaving to higher luminosity 
bins. 



tion models (De Lucia et al., 2006, Benson & Devereux, 2010) 
to study the number density of ETGs of similar PEL and simi- 
lar stellar masses, respectively. In this section we will integrate 
the ETG LFs over two PEL ranges (shown by the dark and light 
shaded areas in Figure [TJ delimited by Msiz = 0) < -20.5, for 
the bright ETG population, and by -20.5 < Mgiz = 0) < -18.5 
for the faint ETG population to obtain their number densities 
ntrightiz) and njamtiz)- The r-band limits can be obtained using 
B - r = 1.15 (Eq.[2). The luminosity Mb{z = 0) = -20.5 corre- 
sponds approximately to a stellar mass of 10" (CDR), a limit 
used by De Lucia et al., in their analysis of semi-analytic galax- 
ies. In this paper, we evolve the r-band z - 0.165 descendant 
luminosities estimated by PIO to z = using Eq. |2] By inte- 
grating over luminosity the method is less sensitive to the exact 
functional form adopted for the LF. Our analysis of mock cata- 
logues in Appendix A2 ensures little influence from systematics 
coming from the LF measurement method or photo-z errors. 

Our measurements of number density evolution from 
MUSYC galaxies will be compared to those already published 
for the datasets COMBO-17, DEEP2 and SXDS, and the De 
Lucia et al. (2006) models of galaxy formation. To study the 
evolution of the number density of ETGs we need z = com- 
parison samples for which we use the z = ETG LF from the 
SDSS measured by B04 for the B-band, and by Be07 for the 
r-band. Figure |2] shows the resulting bright and faint ETG num- 
ber density evolution (top and bottom panels). As can be seen, 
the number density of bright ETGs (top panel) for COMBO-17, 
DEEP2 and SXDS is roughly consistent with that of the present- 
day bright ETGs up to redshifts z ~ 0.7, indicative of little evolu- 
tion since then (CDR). The dashed lines correspond to the evolu- 
tion of model ETG number densities inferred by De Lucia et al. 
(2006) which, for the bright ETGs, lie below the observational 
points for most of the explored redshift range; at low redshift 
where the volume of the samples is small, there is a marginal 
agreement between the De Lucia et al. results and the observa- 
tions (z < 0.35). 

The squares (soUd correspond to the r-band, open to the 
B-band) in this figure show our results for the MUSYC ETG 
LFs, obtained via the PML method. As can be seen, the PML 
measurements in the B-band show a good agreement with all 
the other observational results presented here from DEEP2, 
COMBO-17 and SXDS. The results from both bands, and in par- 
ticular from the r-band, also show a very good agreement with 
the De Lucia et al. (2006) semi-analytic model results for the 
evolution of the bright ETGs. The agreement for the different 
bands varies slightly probably due to the slightly redder colours 
of the MUSYC ETGs in comparison to those from SDSS. 

As can be seen, to the level of certainty allowed by the com- 
bination of the different observational estimates, it is difficult to 
ascertain whether the observed evolution of the bright ETG num- 
ber density is consistent with this particular model. Notice, how- 
ever, that results from some of the catalogues shown here taken 
individually, could be used to rule out this particular model; this 
shows the importance of using different pencil-beam directions 
to be able to account for the important cosmic variance in the 
counts. This also affects to some degree the conclusions from 
Hopkins et al. (2008), since the observational results from z < 1 
only provide loose constraints to disentangle between different 
star-formation quenching mechanisms. 

The evolution of the faint ETG number density shows consis- 
tent results between the different catalogues, and are slightly be- 
low the semi-analytic estimates, even at redshifts below z = 0.65 
where the recovery of the number densities in the mocks was 
succesful (See Appendix A2). However, at least for MUSYC 
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Fig. 2. Ratios of ETG number densities at a given redshift with 
respect to a z ~ SDSS sample. The horizontal solid lines show 
the unit ratio. The shaded areas show the cosmic variance as 
calculated in Appendix A2. The squares show the results from 
MUSYC for the r-band (solid symbols) and the B-band (open 
symbols, slightly displaced on the x-axis) from PML estimates 
of the LF; the errorbars for the results in the two bands (r and B) 
are similar, but are shown only for the r-band to improve clarity. 
The crosses show the results from the r-band LF calculated us- 
ing the photo-z based maximum likelihood method. The circles 
correspond to results from the COMBO- 17 ETG LF estimate, 
triangles to DEEP2, and stars to the SXDS. Top panel: samples 
of bright ETG galaxies. Lower panels: faint ETGs. Comparing 
the shaded area and symbols, a marginal agreement between the 
evolution in models and observations can be noticed. 



galaxies, the expected cosmic variance (See Appendix A2) is 
able to account for some of these differences. 



4. Clustering selection of descendant samples 

In this section we present a new method to select samples 
of galaxies at different redshifts that follow a progenitor-to- 
descendant relationship in a statistical way. Such samples can 
be used to constrain the role of mergers in the evolution of ETGs 
by comparing their number densities to that of their inferred de- 
scendants. This method can also be applied to samples of galax- 
ies with available stellar masses, instead of PELs. 

Instead of adopting a stellar mass selection, the clustering- 
selected descendants are obtained by taking a sample of haloes 
from the simulation with the dark-matter masses correspond- 
ing to the hosts of ETG galaxies at a given redshift obtained 
via clustering measurements (PIG). Then these haloes are fol- 
lowed to z = using merger trees extracted from the numerical 
simulation, and their clustering is compared to that of ETGs in 
the SDSS (z=0, from Swanson et al., 2008). The SDSS ETGs 
with this clustering are then selected as the descendants of the 
MUSYC ETGs. PIO used this method to connect different ETG 
populations at 0.25 < z < 1.4 to ETG galaxies at z < 0.25. This 
study shows that in a ACDM cosmology it is not correct to con- 
sider ETGs with equivalent passively evolved z = luminosities 
at different redshifts as descendants/progenitors of one another 
since the low redshift galaxies with the same PEL/stellar mass 
show a lower clustering than is expected for the descendants. 



Table 1. Properties of clustering-selected z = descendants of 
Mr(z = 0) < -19.7 (stellar masses above ~ lO'^h^^Mo) progen- 
itors at different redshifts. 



Progenitor 
Redshift 



Descendant 

M,{z = 0) 



Descendant 

logio(M,/h-'Me) 



0.175 
0.325 
0.5 
0.7 
1.0 



-17.5; 
-18.8;; 
-19.4; 
-20.2 
-20.8 



Q 4+0-4 

10 r''- 

10 4+" - 

-92 

10.9!; 
11.2!- 



PIO find that MUSYC ETGs with passively evolved z = lu- 
minosities Mriz = 0) < -19.7, evolve from typical luminosities, 
expressed in units of the ETG L* of L/L*{z) - 0.7 at z = 1.15 
into galaxies with typical L/L*{0) ~ 2.1 - 5.2 at z = 0, and from 
L/L*iz) ^ 1 galaxies at z = 0.35 to L/L*{Q) ~ 0.35-0.80 Table 
[T] shows the typical descendant z = B-band magnitude and 
stellar mass for our samples of progenitor galaxies. Notice that 
as progenitors are selected using the same PEL (or stellar mass) 
cut, they are not expected to lie on the same evolutionary line; 
this is confirmed by the different typical descendant properties 
shown in the table. 

We calculate the number density of descendants using the 
B04 and Be07 SDSS ETG LFs for galaxies with median lu- 
minosities corresponding to the descendants of a given sam- 
ple of MUSYC ETGs at redshift z. The number density of 
progenitors is calculated at redshift z, using a lower limit in 
Mriz = 0) = -19.7, or equivalently, Mb{z = 0) = -18.55, 
consistent with the lower limits used in PIO; according to the 
relation between z = B-band luminosity and stellar mass by di 
Serego Alighieri et al. (2005, Figure 13 in their paper), this lu- 
minosity cut effectively selects ETG galaxies with stellar masses 
Ms > lO'^h^'M©. The ratios between the number density of pro- 
genitors and descendants are shown as squares in the upper panel 
of Figure |3] (open symbols for the B-band, filled symbols for the 
r-band). The errorbars in this panel correspond to the uncertain- 
ties in the descendant luminosity, extracted from PIO. As can be 
seen, the ratio is significantly higher than unity at z > 0.6 in both 
bands, indicating the need for mergers between ETGs in order 
to diminish their number density towards z - OQ Lower redshift 
samples show number densities similar to their expected descen- 
dants. As can be seen, the photo-z based maximum likelihood 
method (crosses, shown only for the r-band for clarity) provides 
results in agreement with those from the PML based LFs. 

Taking advantage of the measured LFs for the high-redshift 
ETG samples, we calculate the ratios between the luminosity 
densities of the high-z ETGs and that of their z = likely de- 
scendants, using the PML LF measurements. This is shown in 
the lower panel of Figure|3]as squares (solid symbols correspond 
to results in the r-band, the open symbols to the B-band). As can 
be seen, regardless of the photometric band the data shows that 
as the redshift decreases the luminosity density of ETGs is con- 
sistent with the unit value (bearing in mind that our different high 
redshift ETG samples are not in the same evolutionary line). The 
lowest redshift point is marginally above the unit ratio, but is also 



- The ranges in descendant L/L*(0) correspond to variations in the 
results between the EHDF-S and ECDF-S fields. 

The mock catalogue analysis in Appendix A shows a possible 
overestimate of a factor 1 .5 at this redshift range, which we will bear in 
mind in our conclusions. 
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Fig. 3. Top panel: Ratio between the number density of pro- 
genitor ETGs with stellar masses M, > lO'^h-iMo at a given 
redshift and their expected z = descendants, as inferred from 
the clustering analysis of PIG. The filled squares show the results 
from using the PML r-band LF estimates; open squares show the 
results for the B-band. Errorbars are only shown for the r-band 
result to improve clarity; errors for the B-band show similar am- 
plitudes. The crosses correspond to the results from the photo-z 
based maximum likelihood method to calculate the LF in the r- 
band. The solid line shows the unit ratio and the grey shaded area 
shows the estimated cosmic variance in a 0.25sq. degree light- 
cone survey divided in slices of Az = 0.1. As can be seen, the 
number density of high redshift progenitors is higher than that 
of their z = descendants; this is not the case for progenitors 
at z < 0.6 and their z = descendants. Bottom panel: ratio 
between the luminosity density of progenitor ETGs and that of 
their z = ETG descendants in the r- and B-bands (filled and 
open squares, respectively). Regardless of the progenitor red- 
shift, the luminosity density ratios shown are consistent with the 
unit value. 



influence requires the use of models that include the evolution 
of dark-matter haloes in a fully non-linear way (merger trees) 
including galaxies, we defer this to a forthcoming paper (Padilla 
et al., in preparation). In this paper we present possible sources 
of these systematic effects, and in some cases, qualitative esti- 
mates. 

On the one hand, we have not taken into account the possi- 
bility that some of the stellar content of ETGs may be lost, for 
instance into the intracluster medium (ICM) due to tydal strip- 
ping, before a merger takes place. Conroy, Wechsler & Kravtsov 
(2007) showed that in order to reproduce the evolution of the 
massive end of the stellar mass function since z = 1, the z ~ 
brightest cluster galaxy-cluster mass relation, and the intraclus- 
ter light luminosity distribution, up to 80 percent of the stellar 
content of galaxies in disrupted substructures (i.e. only a frac- 
tion of the satellites) within dark-matter haloes should be lost to 
the intracluster medium; in the next Subsection the analysis of 
the halo model shows some evidence for sinks consistent with 
this estimate. Additionally, the star-formation activity could be 
reignited via new cold gas infall or AGN shutdown, which could 
also remove the galaxy from the ETG sample. On the other hand, 
blue galaxies could become ETGs and contaminate the descen- 
dant population. 

Notice that the approach followed by PIO to obtain a typi- 
cal host DM halo mass does not need information on the stellar 
masses of galaxies, only on their clustering. However it is nec- 
essary to take into account the clustering amplitude vs. stellar 
mass relation since the PIO approach to find descendants re- 
quires a one-to-one relation between PEL or stellar mass and 
host halo mass. In the case of z = 1 MUSYC ETGs this relation 
shows a minimum at Mb{z - 0) - -18.5 (PIO) which indicates 
that our approach can only be applied to galaxies brighter than 
this lower limit. In order to follow the descendants of fainter 
galaxies it would be necessary to adopt, for instance, a subhalo 
abundance matching method. The latter provides a direct con- 
nection between DM subhalo mass (at infall) and stellar mass 
(e.g. Kravtsov et al., 2004), but is subject to uncertainties even 
when clustering and mass function measurements are available 
(Neistein et al., 2001). 



aflected by the largest poisson uncertainties in the LF measure- 
ment due to the small number of galaxies at such low redshifts in 
the ECDF-S, and also by the large effect from cosmic variance 
due to the small volume. 

The results from this analysis show that ETGs at z = 1 are 
likely to descend into z = ETGs undergoing a decrease in 
space-density of a factor 5.5 + 4.0 (combining the results from 
the B and r bands). But, due to the constant luminosity density, 
mergers would provide the required increase in luminosity with- 
out the need for important episodes of star-formation. Notice that 
the amount of mergers derived here (of about four mergers since 
z = 1 for each z = galaxy) is higher than that estimated from 
close pairs of stellar mass-selected samples at different redshifts; 
Robaina et al. (2010) find that galaxies undergo ~ 0.7 merg- 
ers between z = 1 .2 and 0. We suggest that a more consistent 
measurement of merger rates using close pairs would require the 
use of samples in a more likely progenitor-descendant relation- 
ship than that provided by a stellar mass selection. The use of 
clustering-selected samples in this case would include brighter 
low redshift samples which could change the resulting merger 
rates. 

It should be borne in mind that our analysis does not con- 
sider sources or sinks, but since a quantitative analysis of their 



4. 1. Clustering-selected samples and halo model analyses 

Further clues on the evolution of ETGs were found by Brown et 
al. (2008) who analysed the Bootes field. Using the Halo Model 
they point out that central ETG galaxies at z = 0, which occupy 
the bright-end of the LF, have only increased in about 30 per- 
cent their stellar mass since z = 1, whereas the stellar masses of 
satellite galaxies have increased by a factor of ~ 3. 

Applying their universal Halo Occupation Model fits (equa- 
tions 12 and 13 in their paper) to our results for the z = 1 
MUSYC and their descendant z = SDSS ETG populations we 
find that in median, (i) ~ 80 + 5 and ~ 94 + 3 percent of the pro- 
genitor^ and descendants (respectively) are central galaxies, (ii) 
The latter allows us to separate centrals and satellites by sharp 
cutoffs in luminosity as a first approximation, with which we ob- 
tain that centrals increase their luminosity by a factor of 1 ■T'^^'^, 
whereas satellites do so by a factor of 2.5^[ j (consistent with the 
Brown et al. estimates), (iii) and that the number density of cen- 
trals and satellites decreases by factors of ~ 4 + 2 and ~ 10 + 7, 
respectively. Taking into account the relative numbers of cen- 

"* We apply the HOD to the unevolved luminosities of the progenitors 
as its parameters are provided for rest-frame ETG magnitudes by Brown 
et al. 
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Fig. 4. Frequency of stellar mass ratios of ETG galaxies in- 
volved in mergers, for the case of z = 1 progenitors with 
Mriz = 0) < —19.7 (i.e. stellar masses above ~ 
and their z = descendants for a total of 4 mergers per de- 
scendant. The solid line shows the sum of all the merger events, 
dashed corresponds to mergers between central galaxies, dotted 
to mergers between a central and one of its satellites. The shaded 
region corresponds to ratios above the limit for a major merger 
(mass ratio of 1 /3 and higher). 



trals and satellites, this is in agreement with the overall factor of 
5.5 + 4.0 decrease in number density of ETGs. 

The percentage of galaxies which underwent major mergers 
has been predicted by semi-analytic models of galaxy forma- 
tion. Kochfar & Silk (2009) show that between 1% and 2% of 
the massive galaxies (with stellar masses similar to those in our 
samples) undergo a major merger in their last Gyr for z <= 1. 

Our results are consistent with the Kochfar & SiUc estimates. 
The decrease in the descendant MUSYC ETG number density 
can be the product of both minor and major mergers (consid- 
ering the latter those with a mass ratio of 1/3 or higher, e.g. 
Lagos, Cora & Padilla, 2008). Taking into account the shape of 
the progenitor ETG luminosity function at z = 1 and separat- 
ing centrals and satellites by sharp luminosity cuts, we make a 
Monte-Carlo estimate of the ratios between luminosities merg- 
ing galaxies. Our average result of a factor of 5.5 + 4.0 decrease 
in the number density of ETGs points to about four mergers per 
galaxy, where we assume that each one increases the brightness 
of the central galaxy. The ratios between merging galaxies are 
shown in Figure |4] where the shaded region corresponds to ma- 
jor mergers and as can be seen, only a small fraction of the total 
merger events (solid line) fall into this category. We find that 
~ 31% of the central galaxies experience a major merger since 
z = 1, with a decreasing frequency towards z = 0; this translates 
into an average ~ 4% of ETGs undergoing a major merger in 
their last Gyr. We also find that 2/3 of the major mergers involve 
a z = 1 central galaxy that became a satellite and then merged 
with a larger central galaxy. 

This simple calculation ignores possible correlations be- 
tween central and satellite galaxy luminosities and also the evo- 
lution of the ETG LF between z = 1 and 0. Therefore, these 



results should only be taken as a rough estimate of the frequency 
of major mergers. A more thorough analysis would require the 
use of conditional luminosity functions of satellite luminosities 
in terms of those of their central galaxies, on the one hand, and 
more precise measurements of LF and clustering of ETGs in this 
range of redshift, on the other An alternative approach is to take 
into account merger trees of dark-matter haloes and to follow the 
evolution of galaxies within them. We will perform this anal- 
ysis in a forthcoming paper where the dark-matter haloes will 
be selected so as to mimic the high-z MUSYC ETG population 
(Padilla et al., in preparation). 

The Halo Model analysis also allows us to infer the total lu- 
minosity in progenitors (adding centrals and satellites together) 
which exceeds that of descendants by a factor of ~ 4'^^^ which 
points to a possible sink in the luminosity of the ETG population 
of similar amplitude to the expected effect of the ICM pointed 
out by Conroy, Wechsler & Kravtsov (2007). Additionally, the 
hint that the decrease in the number density of satellites is larger 
than that of centrals would indicate that a majority of the merg- 
ers would occur with central galaxies of smaller mass, progen- 
itor dark-matter haloes before joining the final, more massive 
halclfl consistent with independent indications pointing to the 
same conclusion (Gonzalez & Padilla, 2009, Porter et al., 2008). 

5. Conclusions 

In this paper we presented a new method to select ETG samples 
that follow a progenitor-to-descendant relationship which uses 
clustering and luminosity or stellar mass function information. 
We used this method to define clustering-selected ETG sam- 
ples from the MUSYC ECDF-S field, which compared to their 
clustering-selected descendants in the SDSS allowed us to study 
the importance of mergers in their evolution. 

Our main results and conclusions can be summarized as fol- 
lows. 

(i) From the analysis of the number density of samples of 
MUSYC ETGs with equal luminosities, passively evolved 
to z = (a procedure that effectively compares samples 
of galaxies with similar stellar masses), we find an evolu- 
tion of the abundance of bright ETGs characterised by a 
steady increase even down to z ~ 0.3, consistent within 
the errorbars with results from COMBO-17, DEEP2, and 
the De Lucia et al. (2006) measurement using semi-analytic 
galaxies. COMBO-17 and DEEP2 show slightly less evo- 
lution in the ETG number density and a lesser agreement 
with the model, showing the importance of numerous pencil- 
beam like catalogues to pin down cosmic variance. We point 
out that, in principle, this selection does not produce sam- 
ples obeying a progenitor-to-descendant relationship (as as- 
sumed in several works, e.g. CDR, Robaina et al., 2010) even 
though it can still be used for meaningful comparisons be- 
tween models and observations. 

(ii) We implemented the clustering selection method with the 
aim of finding clues on the role of mergers in the evolution 
of ETGs. We define progenitor ETG samples from MUSYC 
with median redshifts between 0.2 and 1. Their z - 0.1 de- 
scendant galaxies are selected so as to have the same cluster- 
ing amplitude as the dark-matter haloes that host the high- 



' The reason for this is that if there were no galaxy mergers in the 
progenitor haloes, most of the z = 1 satellites would survive until today 
due to the fact that the dynamical friction timescales are much larger for 
lower satellite-to-central mass ratios (see for instance Boylan-Kolchin 
et al., 2008). 
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z ETG samples (obtained from clustering measurements by 
PIO) followed down to z = in numerical simulations. The 
analysis of these samples allowed us to find some evidence 
for dry mergers of ETG galaxies between redshifts z ~ 0.8 
and z ~ 0.2, since their number density declines by a factor 
~ 5.5 +4.0 (Figure|3] top panel, the ranges show the variation 
in the results from using the B- and r-bands). Possible sys- 
tematics in the LF measurement obtained in our analysis of 
mock catalogs could lower this factor to ~ 3.8 + 2.6. Given 
that the luminosity density of progenitors and descendants 
are consistent with one another (Figure |3] bottom panel), 
our results point towards dry mergers without important star- 
formation episodes in ETGs, at least down to the precision 
allowed by the MUSYC dataset. These two results are con- 
sistent in the sense that at z = there are not enough ETGs 
to account for the progenitor population via simple passive 
evolution; the merging scenario helps solve this problem, by 
decreasing the number density of fainter galaxies, which af- 
ter merging provide the observed ETG luminosity and num- 
ber density at z = 0. 
(iii) Combining our results for the descendants of z = 1 MUSYC 
ETGs with halo occupation fits to LF and clustering mea- 
surements by Brown et al. (2008), we find that ~ 4% of the 
descendants of z = 1 ETG galaxies with stellar masses above 
~ 10'"h"'Mo undergo a major merger in their last Gyr, in 
agreement with recent results from semi-analytic models by 
Kochfar & Silk (2009). We also find that two thirds of the 
major mergers involve an ETG which was central galaxy at 
z = 1, and some indications that most of the ETG mergers 
would need to take place in groups before these fall into large 
clusters (in agreement with Porter et al., 2008, and Gonzalez 
& Padilla, 2009). 

In the analysis of (ii) and (iii) we have assumed that the 
Z = ETG sample with the same clustering as the descendants of 
the haloes hosting high redshift ETGs does not include galaxies 
which turned red and dead in between. These would contribute to 
the number density of ETGs at z = 0, which would then be over- 
estimating that of the true descendants of the individual ETG 
progenitors. On the other hand, our approach has not considered 
sinks for the ETG population, which could be provided by tidal 
effects in the intracluster medium or by processes that reignite 
the star-formation activity in galaxies, which in turn would over- 
estimate the number density of the progenitors. 

Upcoming surveys such as those planned for the Large 
Synoptic Survey Telescope (Abell et al.. The LSST Science 
Book, 2009), will allow much better statistics by increasing the 
solid angle with respect to the currently available deep photo- 
metric surveys such as MUSYC and COMB017. Furthermore, 
with a higher signal it could also be possible to produce measure- 
ments of the expected merger rates and star formation history of 
combined samples of early and late type galaxies, and via com- 
parisons with models of galaxy formation to help improve our 
understanding of the evolution of galaxies from high redshifts to 
the present day. 
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Appendix A: Mock catalogue analysis 

We use the ECDF-S mock galaxy catalogues from C09 (two in 
total) to analyse the effects of distance uncertainties on the LF. 
The mocks are extracted from a ACDM numerical simulation 
populated with GALFORM (version corresponding to Baugh et 
al., 2005) semi-analytic galaxies. The simulation contains 10^ 
DM particles in a periodic box of lOOOh 'Mpc a side. The cos- 
mological model adopted in this simulation is characterised by a 
matter density parameter Q,„ = 0.25, a vacuum density parame- 
ter Qa - 0.75, a Hubble constant, H - /zlOOkms 'Mpc ', with 
h = 0.7, and a primordial power spectrum slope, n , - 0.97. The 
present day amplitude of fluctuations in spheres of Sh^'Mpc is 
set to (Tg = 0.8. This particular cosmology is in line with recent 
cosmic microwave background anisotropy and large scale struc- 
ture measurements (WMAP team, Dunkley et al., 2009, Spergel 
et al. 2007, see also Sanchez et al., 2006). We adopt this cosmol- 
ogy for all the calculations performed throughout this paper. 

In order to make the mocks as similar to the observations 
as possible, C09 used the information on stellar mass and bulge 
fractions to assign a SED drawn from a smooth continuum of 
possible templates (from the same set used to calculate photo- 
zs) to each galaxy. This SED is used to calculate fluxes in the 
MUSYC bands. Using these photometric data, C09 repeated the 
photo-z measurement procedure on the mock galaxies, which 
are affected by similar Zphot-Zspec eiTor distributions as the actual 
MUSYC data, including outliers and catastrophic eiTors. 

The two mock catalogues match the angular selection func- 
tion of MUSYC fields but include no evolution since they were 
constructed using only the z = simulation output. In order to 
mimic an active evolution of the mock ETG galaxy luminosities 
with redshift, we apply a linear magnitude brightening of 0.6 
magnitudes between z = 1 and z - 0. This choice produces an 
evolution of the number density of bright ETGs consistent with 
the reported evolution of ETGs in the De Lucia et al. (2006) 
model (also consistent with the evolution of GALFORM ellipti- 
cals found by Benson & Devereux, 2010). We apply this bright- 
ening of the LF to our mock results from this point on. 

In this appendix we use six extimates of the LF. Those ob- 
tained using the LF of ETG galaxies in the full simulation box, 
the PML and classical photo-z based maximum likelihood esti- 
mates, and 1 /Vmax estimates of the LF from the mocks using the 
true galaxy redshifts (with no errors), the photometric red shifts, 
and the latter with a fuither convolution with the photo-z errors. 

Figure lA.ll shows the LF of the combined set of galaxies 
from both mock catalogues; this is done in order to minimise 
random errors and focus our analysis on possible systematic ef- 
fects. From left to right, the panels show different redshift slices 
selected by placing limits on photo-zs (the average redshifts are 
indicated in each panel). The triangles represent the true, under- 
lying luminosity functions in the simulation, and circles show 
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Fig. A.l. Effect of redshift measurement errors on the luminos- 
ity function as found using mock MUSYC catalogues, analysed 
in three redshift bins (increasing from left to right). In all the 
panels, the triangles show the tiTie, evolving r-band ETG LF of 
the full simulation box. Top panels: LFs calculated using the 
1/Vmnx weighting method. Filled circles show the results from 
using the true redshifts in the mock catalogues, the solid lines 
indicate the results when adopting photo-zs, the dashed lines 
shows the effect of a further convolution with the photo-z er- 
rors. Bottom panels: the filled circles are as in the top panels. 
The solid lines coiTespond to the PML estimates of the LF from 
the mock, the dashed lines to the photo-z based maximum like- 
lihood estimator (C09). Shaded areas and vertical dashed lines 
are as in Figure [1] The vertical dotted line shows the magnitude 
above which the recovery of the LF using the true redshift starts 
to degrade. EiTorbars show poisson uncertainties. 



-21.5 




-22.5 - 



Fig. A.2. Recovered values of M* from the LFs in the top panels 
of Figure lA.ll as a function of redshift. The results from fitting 
the true LF are shown as triangles, from the true redshifts in the 
mock as circles, from the photo-zs as solid lines, and from the 
double convolution as dashed lines. Errorbars show the 1 - cr 
confidence ranges obtained from the fits. Small displacements in 
redshift were applied to improve clarity. 
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the l/V,„ax estimate of the LFs from the mocks using the true 
galaxy redshifts. 

In the top panels the lines correspond to LFs calculated us- 
ing the 1 / Vrnfl v weighting with photo-zs. The solid line shows the 
resulting LP from assuming the photo-z provides the correct dis- 
tance to the galaxy; this approach produces an excess of bright 
galaxies (with respect to the result from the true redshifts) which 
is slightly more important at lower redshifts. If in addition to the 
photo-z, a convolution of the distance with the photo-z error is 
applied to the data, the bright-end of the LF is once more sys- 
tematically enhanced (dashed lines) and can translate into a shift 
of the LF towards brighter luminosities of up to one magnitude. 
The algorithms applied to the SXDS and COMBO- 17 data are 
comparable to the latter approach, and therefore could produce 
results biased towards a high space density of very bright ob- 
jects. Notice, however, that most studies of the luminous ETG 
population integrate the LF over a wide range of luminosities 
marked by the dark or light-grey shaded regions, and therefore 
the effect on such statistics could be only minor; in particular, 
the number density of galaxies brighter than M* are underesti- 
mated in the two photo-z based 1 / Vmax LF measurements, by 
only ~ 0.2dex (see also the following Section). 

The solid and dashed lines in the bottom panels, show the 
results from the PML and photo-z based maximum likelihood 
methods (respectively). As can be seen, in comparison to the 
1 /Vmax estimates from photo-zs, the photo-z based maximum 
likelihood method shows a much milder tendency of shifting 
the bright-end toward brighter magnitudes (with respect to the 
PML method and 1/V,„ax tme-redshift estimates which pro- 
vide the best and second-best matches to the true LF, respec- 
tively). This brightening is almost completely counter-balanced 
by a slight decrement in the space-density of faint galaxies. The 
PML method provides a good match to the mock result ob- 
tained using the true galaxy redshifts, and is marginally better 
than the photo-z based maximum likelihood result. Regardless 
of the measurement method, the recovery of the faint-end of the 
LF starts to diverge from the true LF at luminosities 0.3 to 0.6 
magnitudes brighter than the completeness limit (shown by the 
dotted and dashed lines, respectively). We take this into account 
when studying the space density of faint ETGs in observational 
catalogues. In particular, for the analysis of Section 4, the space 
density of z > 0.6 ETGs with Meiz = 0) < -18.55 could be 
overestimated by 50%. 

A.1. Recovery of the characteristic luminosity, M* 

Figure lAT2] shows the recovered values of the characteristic lumi- 
nosity M* from the LFs of the full simulation, the mocks using 
true redshifts, using photo-zs, and from the double convolution 
with errors. The evolution of 0.6 magnitudes per unit redshift is 
recovered in all cases and the LF from the true-redshifts is con- 
sistent with the true underlying LF. The brightening of the LFs 
affects the recovered values of M* for the case of photo-zs and 
double convolution of ^ 0.6 and 1 magnitudes at z = 1.1, respec- 
tively. This effect becomes slightly less severe at lower redshifts. 

A.2. Number density evolution with mock catalogues 

We measure the evolution of the ETG number densities using the 
LFs from the mock catalogues to analyse the systematic effects 
arising from the use of photo-z based LF measurements. 

Figure IA.3I shows the ratios between the number density of 
mock ETG galaxies at a given redshift with respect to the z = 
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Fig.A.3. Ratios of ETG number densities at a given redshift 
with respect to the z ~ true luminosity function in the numer- 
ical simulation. The horizontal solid line shows the unit ratio. 
The solid squares show the results from the mocks for the r-band 
from PML estimates of the LF; the errorbars are shown only for 
this case to improve clarity (and have been scaled to an area of 
0.25sq. degrees). The crosses show the results from the photo- 
z maximum likelihood measurement. The open symbols show 
the results from the 1 /Vmax weighted estimates of the mock LFs, 
obtained using the true redshifts (open squares), the photometric 
redshifts (open circles), and the photo-zs convolved with their 
uncertainties (open triangles). Small displacements along the x- 
axis have been added to improve clarity. The dashed lines cor- 
respond to the evolution inferred for the De Lucia et al. (2006) 
semi-analytic model (coincident with the evolution enforced in 
our model galaxies; see the text). The grey shaded area shows 
the estimated cosmic variance in a 0.25sq. degree light-cone sur- 
vey divided in slices of Az = 0.1. Top panel: samples of bright 
ETG galaxies. Lower panels: faint ETGs; the dotted vertical line 
shows the redshift limit for accurate faint ratios. 

ETG LF from the full simulation box, mimicking the analysis 
of Section 3. The dashed lines embedded in the shaded areas 
show the underlying evolution of this ratio and the expected cos- 
mic variance for the adopted cosmology, for spherical volumes 
equivalent to Az = 0. 1 thick redshift slices of light-cones cover- 
ing a solid angle of 0.25sq. degrees (appropriate for the ECDF- 
S). In this calculation we use the Smith et al. (2003) approxima- 
tion to the non-linear power spectrum of density fluctuations. 

As can be seen in the figure, both the PML and photo-z max- 
imum likelihood methods provide a good match to the underly- 
ing evolution, with a slight overestimation at redshifts z > 0.7 
which is more severe for the faint population, as expected from 
the analysis of the previous section. Errorbars represent the un- 
certainties in a field of 0.25sq. degree area. The other symbols 
show the l/Vmiii estimates of the LF obtained from the true 
galaxy redshifts (the best fit to the evolution in the full simula- 
tion box), photometric redshifts, and from photo-zs which have 
been additionally convolved with the photo-z estimate errors. As 
can be seen, in the latter two cases, the bright ETG number den- 
sity systematically higher although still in agreement with the 
underlying evolution. 

The figure also shows that for faint galaxies the systematic 
effect from redshift errors is negligible compared to flux incom- 
pleteness effects which become important at these luminosities. 
Even if the amplitudes of the photo-z errors were considerably 
larger for fainter objects this would still affect the faint ratios 
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relatively little unless the faint-end of the LF was very steep 
(i.e. a « -1). Even though, as was shown in Figure IaTI 
the flux incompleteness limit lies within the limits of the faint 
population for the highest redshift bin, there are important dis- 
crepancies with the true faint ETG number densities for z > 0.7. 
Furthermore, the parametric fits to the LF overestimate the faint 
ETG densities, whereas the 1 /Vmax method underestimates them, 
indicating an important difference between the two methods. 



